Analysis of the Kondo effect in ferromagnetic atomic-sized contacts 
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Atomic contacts made of ferromagnetic metals present zero-bias anomalies in the differential 
conductance due to the Kondo effect. These systems provide a unique opportunity to perform a 
statistical analysis of the Kondo parameters in nanostructures since a large number of contacts 
can be easily fabricated using break-junction techniques. The details of the atomic structure differ 
from one contact to another so that a large number of different configurations can be statistically 
analyzed. Here we present such a statistical analysis of the Kondo effect in atomic contacts made 
from the ferromagnetic transition metals Ni, Co and Fe. Our analysis shows clear differences between 
materials that can be understood by fundamental theoretical considerations. This combination of 
experiments and theory allows us to extract information about the origin and nature of the Kondo 
effect in these systems and to explore the influence of geometry and valence in the Kondo screening 
of atomic-sized nanostructures. 



I. INTRODUCTION 

The Kondo effect is one of the most intriguing phenom- 
ena arising from electronic correlations which was first 
observed over 80 years ago as a then unexpected increase 
of the resistance of gold wires at very low temperatures^ 
This phenomenon was successfully explained by Kondo 
30 years later in his seminal work^ as due to scattering 
of conduction electrons off magnetic impurities present in 
the Au samples, thereby screening its magnetic moment. 
More generally, whenever a local magnetic moment is 
coupled to a sea of conduction electrons, the Kondo ef- 
fect can arise at low enough temperatures with important 
consequences for the electronic and magnetic properties 
of the system. 

Also in the case of mesoscopic devices, the Kondo ef- 
fect strongly alters the electronic structure and therefore 
has dramatic consequences on the transport characteris- 
tics of the system. One of the simplest mesoscopic de- 
vices showing the Kondo effect is the case of a quantum 
dot connected in series to two metallic leads; 3 '^ When 
an electronic level of a quantum dot is well below the 
Fermi Energy of the metallic leads and the Coulomb re- 
pulsion is strong enough to prevent double occupation, 
the quantum dot behaves as a magnetic impurity. In this 
situation conduction through the quantum dot is usually 
strongly suppressed due to the Coulomb blockade^. How- 
ever, at low temperatures the Kondo effect restores the 
conductance due to the apparition of a sharp resonance 
-the so-called Kondo resonance- in the spectral function 
of the quantum dot right at the Fermi energy of the elec- 
trodes. 

The essence of the Kondo effect is the formation of 
a total spin-singlet state between the impurity electrons 
and the conduction electrons near the Fermi level^ be- 
low a certain critical temperature characteristic of the 
system, the Kondo temperature. The formation of this 



Kondo singlet state gives rise to the effective screening of 
the magnetic moment of the impurity, and leads to the 
formation of a sharp resonance in the spectral density of 
the impurity electrons right at the Fermi level. This is 
the affore mentioned Kondo resonance, sometimes also 
called Abrikosov-Suhl resonance^—. In the case of mag- 
netic impurities in metallic host materials, the formation 
of the Kondo resonance in the spectral density of the 
impurity leads to additional scattering of the conduction 
electrons, resulting in the increase of the resistance of the 
metal at low temperatures (for a review of the Kondo ef- 
fect in bulk metals with magnetic impurities see e.g. Ref. 

Other systems where evidence for the Kondo effect has 
been found by the manifestation of a Kondo resonance 
either in the spectral density or the conductance charac- 
teristics, include point contact s 10 ^ 1 , different molecules 
containing magnetic atoms on surface o 12 ' 13 , fullerenesi^, 
carbon nanotubes contacted by metallic electrode o 15 i 16 
and magnetic atoms on surfaces studied by Scanning 
Tunneling Microscopy (STM)iS - ^ 

In the case of magnetic adatoms on metal surfaces 
studied by STM, interference of different conduction 
channels through the atom (one of them bearing the 
Kondo resonance), gives rise to Fano lineshapes 2 ^ in the 
low-bias conductance characteristics, similar to the case 
of a quantum dot coupled laterally to a wirei 27 ' 28 By fit- 
ting those lineshapes to the Fano model one can obtain 
different parameters that describe the characteristics of 
the Kondo screening in the system, i.e. the width, po- 
sition and amplitude of the Kondo peak (for a review 
see e.g. Ref. l20h . Of special relevance to our work 
are the STM experiments performed in the high conduc- 
tance regime when the tip is brought into contact with 
the adatom 2 ^— . 

In a recent work we reported the observation of Kondo- 
Fano lineshapes in the conductance characteristics of 
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atomic contacts made from ferromagnetic materials^. In 
contrast to 

STM experiments where the contacted adatom can be 
imaged and the geometry of the system formed by the 
adatom and the surface can be completely characterized, 
break junction experiments do not allow to control the 
geometry of the system to the same extend (although nor- 
mally give rise to some geometrical repetition^). On the 
other hand, using break junctions it is possible to form 
and study a large number of different configurations^ 

The observation of the Kondo effect in ferromagnetic 
atomic contacts was highly unexpected for two reasons: 
First, given the chemical homogeneity of the atomic con- 
tacts, a devision into magnetic impurity and a spin- 
degenerate Fermi sea, as required for the appearance of 
the Kondo effect, is not obvious at all. Second, electrodes 
and the contact atoms are made from a ferromagnetic 
material. However, Kondo effect and ferromagnetism are 
generally competing phenomena: For example, a strong 
enough magnetic field corresponding to a Zeeman energy 
above the binding energy of the Kondo singlet (about the 
Kondo temperature) will break up the singlet and thus 
the Kondo effect 4 ^ while, on the other hand, magnetic 
fields on the order of the Kondo temperature or below 
will lead to a s plit ting of the Kondo resonance as shown 
e.g. in Refs. HlH and@. 

Likewise the Kondo resonance splits when a quantum 
dot is connected to ferromagnetic electrodes in the case of 
parallel alignment of the two electrodes' magnetic polar- 
izations while for antiparallel alignment of the electrodes 
a normal Kondo effect is obtained^—. Therefore one 
would expect that the coupling of the contact region to 
the strongly ferromagnetic bulk electrodes should either 
eliminate the Kondo effect completely or at least split the 
Kondo resonance unless the antifcrromagnetic coupling 
between magnetic impurity and conduction electrons is 
strong enough compared to any other interactions. 

Here as in our previous works 8 ., we propose that the 
Kondo effect in ferromagnetic atomic contacts origi- 
nates from individual d-levels of the undercoordinatcd 
tip atoms of the nanocontact, and the Kondo screening 
is due to the delocalized sp-electrons which are basically 
spin-unpolarized. Hence the impurity or quantum dot 
bearing the spin consists of one or several (i-levels of an 
individual contact atom. 

In our model, depending on the material the spin of 
the tip atom can be localized in different (i-levels. Due 
to the low symmetry of the contact atoms the individual 
(i-levels couple differently to the sp-conduction electron 
bath, resulting in different Kondo screenings of the spin. 
This scenario may explain the different Kondo behaviour 
observed for contacts made from Fe, Co and Ni.— Con- 
tacts made of the same material may also present slightly 
different atomic configurations which will also influence 
their Kondo properties. In our experiments, from the 
fitting of the differential conductance curves performed 
on atomic-sized contacts of ferromagnetic materials to a 
Fano expression, we extract the relevant parameters that 



characterize the Kondo effect on each contact. 

In this work, the validity of the Fano-Kondo model 
for the contact regime and more specifically for the case 
of atomic contacts is revised and discussed. We present 
a more exhaustive analysis than in Ref. HI for all the 
parameters of the Fano-Kondo model in our system and 
extract new information from the data, for instance, the 
distributions for the occupation of the d-level for Fe, Co 
and Ni. This further analysis confirms the marked dif- 
ferences of the case of Ni with respect to Fe and Co. We 
propose here a new interpretation of the data: the nature 
of the Kondo screening is different for Ni contacts than 
for Fe and Co. From the combined statistical analysis of 
the Kondo parameters and theoretical fundamental con- 
siderations, we can deduce the influence of valence and 
environment on the Kondo screening in nanostructures. 

This paper is organized as follows: The first two sec- 
tions are introductory. First, we review the transport 
properties of atomic contacts in Sec. [Til Then in Sec. IIII| 
we present the basic elements of the theory of the Kondo 
effect in the framework of the Anderson impurity model 
and justify the validity of the Fano-Kondo model for the 
case of atomic contacts of transition metals. Sees. IIVI 
and IVl are devoted to the experimental methods and the 
presentation and discussion of the experimental results 
respectively. We finally discuss the results in the frame- 
work of the theory presented in Sec. I VII and summarize 
our main conclusions in Sec. IVIII 



II. TRANSPORT IN ATOMIC CONTACTS 

The system under study in this work are atomic 
contacts^, i.e. a contact between two metallic leads me- 
diated by an atom forming a metallic bond with both 
leads. More specifically, we focus here on the case of ho- 
mogeneous contacts, those where all the atoms forming 
the structure, leads and contact, are of the same element. 
The formation of an atomic contact can be identified from 
electronic transport measurements. When two pieces 
of the same pure metal are brought into contact, the 
conductance plotted against the inter-electrode distance 
shows a plateau when the atomic contact is formed, close 
to the quantum of conductance Go = 2e 2 /h. The exact 
value of this quantity will depend on the material and 
the atomic configuration of the contact. In general, the 
atomic orbitals will define a number of eigen- channels, 
each of these with a transmission probability Ti which 
will be reflected in the total conductance through the 
Landauer formula G = Gq T;. The number of chan- 
nels will be related to the valence of the meta l 48 ' 49 and 
will normally include a high transmitive s-type channel 
and several other with lower transmission, resulting in a 
final conductance in the range of 0.7-2.5 Go- Experimen- 
tally the evolution of the conductance can be recorded 
over the formation or breaking process of the contact 
showing in most of the cases^ plateaus, not only for the 
one-atom contact but for every atomic rearrangement of 
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the wire while pulling. 
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FIG. 1. (a) Typical traces of conductance for breaking con- 
tacts of Co, Fe and Ni. (b) Histograms of conductance are 
constructed from thousands of these breaking traces of con- 
tacts fabricated by STM at 4.2 K for Ni, Fe or Co. The 
histograms for the three materials show a clear first peak at 
values above 2e 2 fh. 



the d-channels are spin-polarized while the s-channel ba- 
sically is unpolarized. 

III. BASIC ELEMENTS OF THE THEORY 

As discussed in the introduction, in the experiments 
the presence of the Kondo effect is reflected in the con- 
ductance of a one-atom contact as a zero-bias anomaly. 
The details of this anomaly can be related to the char- 
acteristics of the Kondo effect in the framework of the 
Anderson impurity model (AIM) £1 The purpose of this 
section is to provide the basic elements for a theoreti- 
cal description of the Kondo effect in the framework of 
the single-level Anderson impurity model (1AIM) and the 
two- level Anderson impurity model (2 AIM). This will al- 
low us to properly analyse our experiments for different 
materials. Additionally, we develop a simple microscopic 
model in order to understand the occurrence of different 
Fano lineshapes within the same material in terms of the 
variation of microscopic interactions due to variation in 
the atomic structure in the contact region. 

A. Model of a nanocontact 



The one-atom plateau of conductance for the case of 
gold is near the quantum of conductance and, as shown 
in Fig. QJ,, higher for the case of the 3d-transition metals 
Fe, Co and Ni. From many of these traces we can build 
histograms which help us to identify the transport prop- 
erties of the most probable atomic configurations of the 
contacts. 

Different author o 51 ' 52 agree that under cryogenic con- 
ditions and considering a high enough amount of data 
(thousands of conductance traces), the conductance of 
monatomic contacts made from Fe, Co or Ni takes a value 
higher than the quantum of conductance Go = 2e 2 /h, as 
expected in general for transition metals. 40 Histograms 
of conductance of Fe, Co and Ni constructed from thou- 
sands of breaking traces at 4.2 K are shown in Fig. [TJd. 
The most probable values for the conductance of the 
monatomic contacts are 1.2 Go, 1.6 Go and 2Go for Co, 
Ni and Fe, respectively. Conductance values between 
1 Go and 2 Go for ferromagnetic nanocontacts are in 
overall agreement with theoretical calculations i^ 3 - - — The 
conduction of these atomic contacts is the sum of the 
contributions of different channels, where the s-channcl 
is expected to be open and practically degenerate in spin 
and thus to have an associated conductance of nearly 
2e 2 /h. In the ferromagnetic metals Fe, Co and Ni another 
contribution to the overall conductance comes from the 
transmission of electrons via the five 3d-orbitals. How- 
ever, due to the directionality and stronger localization 
of the d-orbitals, the electrons in the d-channels are eas- 
ily scattered. Hence the transmission in d-channels is 
usually far from perfect as has been shown by ab initio 
calculations 53 . In this work it was also shown, that only 



We assume that the Kondo effect takes place in the un- 
der coordinated tip atoms of the nanocontact. This as- 
sumption seems reasonable considering that the Kondo 
effect is not observed in bulk samples of the ferromag- 
netic materials studied here. Hence the Kondo effect 
must be related to the atomic-size constriction of the 
nanocontact. Furthermore, the Coulomb interaction and 
the localization of the electrons within the atomic-size 
constriction, and especially at the tip atoms should be 
enhanced as compared to bulk. 

Typically, the nanocontacts of the ferromagnetic met- 
als Fe, Co and Ni form dimers in the last step before 
breakin g 58 ' 59 as illustrated by the cartoon in Fig. HJa). 
In his case each tip atom will be more strongly coupled to 
one electrode than to the other. As was said before in Sec. 
III1 the main conduction channel through the tip atoms is 
the spin-degenerate s-channel with nearly perfect trans- 
mission. The d-channels on the other hand are strongly 
spin-polarized and their transmission is weaker or even 
completely blocked due to scattering by the geometry, 53 
As we have shown in our previous work 3 ^ by means of 
ab initio calculations, due to disorder and low coordi- 
nation in the contact region one of the ci-levels of the 
tip atom couples only very weakly to the (i-levels of the 
neighbouring atoms, and instead only couples to the ba- 
sically spin-degenerate s-channel. Now this is the situa- 
tion where the Kondo effect can take place and which is 
described by the Anderson impurity model: A strongly 
interacting d-level couples to a non-interacting sea of con- 
duction electrons. Hence we can identify this d-level with 
the "quatum dot" or impurity in our experiments. This 
situation is schematically depicted in Fig. fJfb) . 
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(a) 




FIG. 2. (a) Cartoon of nanocontact just before breaking. 
The tip atom (red) of the left electrode is only weakly coupled 
to the right part of the nanocontact (and vice versa). This 
leads to the simplified model of the nanocontact as shown 
in (b). In the simplified model the tip atom (red) consists 
of a non-interacting s-level and a strongly interacting d-level 
which couple more strongly to the left electrode than to the 
right electrode (Tt,,ct > Yb,,c)- The hybridization V s d between 
s- and d-level is due to crystal field splitting and is weak 
compared to the other energies. Also the coupling of the d- 
level to the electrodes is much weaker than the coupling of 
the s-level (T B > T d ). 



This scenario is further supported by recent Dynami- 
cal Mean-Field Theory calculations of a Ni nanocontact 
connected to electrodes made from Cu instead of Ni, 
thus neglecting the ferromagnetic coupling to the bulk 
electrodes^ In this situation a Kondo effect emerges in 
one of the d-channels of the tip atoms of the Ni nanocon- 
tact with the Kondo temperature in good agreement with 
the ones measured for Ni nanocontacts. 



B. Kondo effect in the Anderson model 

As we have argued in our previous work on the basis of 
ab initio calculations the combined effect of under coor- 
dinated tip atoms and disorder in the contact region can 
lead to the selection of a single d- level of the tip atom that 
is only weakly coupled to the d-levels of the neighbouring 
atoms. Hence the magnetic coupling to the neighbouring 
atoms is reduced for this level and therefore the Kondo 
effect can arise if this level is half-filled. This situation 
can be fulfilled in the case of Ni nanoconctacts where the 
rf-shell is d 9 and thus has a hole. In this case the single- 
level Anderson impurity model (1AIM) where a single 
strongly interacting d-level is coupled to a bath of non- 
interacting conduction electrons is a good description of 
the situation: 

Uaim = £dh d + U hd-fhdi + 2J £ q c\ a c qa 



(i) 



q.a 



where €d is the energy of the d-level, U is the (effective) 
Coulomb repulsion between two electrons in the d-level, 
t q is the energy dispersion of the bath electrons, and V q 
is the coupling (or hopping) between the bath and the 
d-level. In the case of Fe and Co there is more than 



one hole in the d-shell and hence one should consider a 
multi-level AIM. We will do so in the next subsection. 

At zero temperature the Anderson model is a Fermi 
liquid and hence has a (renormalizcd) quasi-particle res- 
onance (i.e. the Kondo peak in the Kondo regime) near 
the Fermi level. The Green's function of the d-level is 
given by: 



G d (u) = 



ui — ex + iY 



K 



(2) 



where z is the quasi-particle weight (i.e. the renormal- 
ization of the single-particle wave-function due to many- 
body effects), e K is the position of the quasi-particle peak 
with respect to the Fermi level (set to zero for conve- 
nience), and Tk = kTx is half the width of the quasi 
particle resonance which defines the Kondo Temperature. 
Correspondingly, the projected density of state of the d- 
level is a Lorentzian centred at £k and width 2Yk ■ 



Pd{uj) 



-ImGd(w) 



zY K /n 



(3) 



In the case of the single-level AIM, the Kondo temper- 
ature Tk can be easily calculated from the parameters of 
the model as follows^ 



r 



K 



K 



^ T d U ire d (e d +U)/T d U 

2 



(4) 



where Yd is the broadening of the <i-level due to the cou- 
pling to the bath obtained by integrating out the bath 
degrees of Freedom: 



Td = ^ 2 pbath(w = 0) 



(5) 



Here we have assumed an approximately constant bath 
density of states pbath and the coupling V independent 
of q. Note the exponential dependence of the Kondo 
temperature on the interaction U and broadening Td- 
This means that mild changes in the parameters can have 
a huge effect on the Kondo temperature. Also note that 
other definitions of the Kondo temperature may differ by 
a constant prefactor. 

Since the 1AIM is a Fermi liquid at zero temperature^ 
we can exploit further relations of the Fermi liquid the- 
ory. For example, one obtains the following important 
relationship between the impurity-level occupation and 
the Kondo parameters (see e.g. Ch. 5 in Ref. @): 



7T \kT K 



(6) 



Also from the Fermi liquid theory of the Anderson model 
we obtain the following exact relation between the den- 
sity of states at the Fermi level tF = (in general not 
the maximum of the Kondo peak) to the occupation of 
the d-level and the broadening Yd due to the coupling of 
the d-level to the rest of the system: 



Pd(0) 



sin 2 (|rt d ) 
TrY d 



(7) 



Finally the amplitude of the Kondo resonance is 

A K = p d {e K ) = — ^— (8) 



C. Underscreened Kondo effect in the multi-level 
Anderson model 

Co and Fe feature 2 and 3 holes, respectively, in the 3d- 
shell of each atom and hence have an atomic spin 5 > 1/2 
due to Hund's rule coupling. Therefore a description in 
terms of a 1AIM as before is problematic. Neverthe- 
less the experimental results can be fitted very well to a 
1AIM (see below and Sec. IV)) . The explanation might be 
that we are really dealing with a so-called underscreened 
Kondo effect^ (UKE) where only a spin- 1/2 in one of the 
d-levels is screened while the rest remains unscreened. 

Such an UKE behaves in many ways like a normal 
(fully screened) 5=1/2 Kondo effect. For example, 
it is still characterized by a zero-bias anomaly resulting 
from a sharp resonance in the screened impurity level. 
Further support for this hypothesis comes from a paper 
by Perkins et al.£^ They find that for an underscreened 
Kondo lattice, ferromagnetism and Kondo effect can in 
fact coexist. On the other hand, the UKE has a number 
of peculiar consequences such as the formation of a so- 
called singular Fermi liquid state characterized e.g. by 
the divergence of the quasi-particle weight and thermo- 
dynamic quantities such as the specific heat capacit y 64 ' 65 . 
Note that in this sense the UKE is very similar to the so- 
called ferromagnetic Kondo effect 66 where the impurity 
spin couples ferromagnetically with the conduction elec- 
trons. The resulting antiscreening of the magnetic mo- 
ment by the conduction electrons also leads to the forma- 
tion of a singular Fermi liquid state. The ferromagnetic 
Kondo effect has recently been discussed theoretically in 
the context of magnetic impurities in nanocontacts i 67 ' 68 

We consider a two-level Anderson impurity model 
(2AIM). This should model the situation of Co which 
has two holes in the d-shell. As in the case of the simple 
AIM, the system is divided into two subsystems: the con- 
duction electron bath B and the impurity I with the two 
interacting levels. The version of the 2AIM model that is 
relevant for the underscreened Kondo effect is depicted 
schematically in Fig. Eta): Only one of the impurity 
levels couples to the conduction electron bath while the 
coupling of the other level is neglig ible^^ Hence the 
Hamiltonian of the 2AIM is given by: 

H2AIM = Hi + Hr + Ht 

= 22 (eifti + U n lf h lir ) + U' hih 2 - J H §1 • s 2 

i 

+ Y e i c U c i° + Y Vi 'i (4 CT <v + c \o d i°) ( 9 ) 

where di„ (d\ a ) destroys (creates) one electron in impu- 
rity level i with spin a, and hi a = d\ a di a is the occupation 
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(a) (b) 




FIG. 3. (a) Schematic drawing of 2-level Anderson impurity 
model in the underscreened situation. Only level 1 of the 
impurity is coupled to the bath B. (b-d) Schematic illustration 
of a hopping process contributing to the underscreened Kondo 
effect in the 2-level Anderson impurity model. 



number operator for level i and spin a and n, = hif+hii- 
Si = J2aa' dl a T aa >di a i measures the spin in level i. ti are 
the energies of the two impurity levels, U is the Coulomb 
repulsion within the same level i and U' is the Coulomb 
repulsion between electrons in different levels which is 
generally smaller than U, and Jh is the Hund's rule cou- 
pling. The bath B is described as in the case of the 
1AIM, eq. (QJ. Finally, only impurity level 1 is coupled 
to the conduction electron bath with hopping V±, while 
the coupling of impurity level 2 is negligible. 

For the sake of simplicity we also assume that the two 
impurity levels are degenerate: e\ — e 2 — e. Typically, 
the intra-level Coulomb repulsion U is bigger than the 
inter-level Coulomb repulsion U' by an amount of the 
order of the Hund's rule coupling: U ~ U' + Jh- Assum- 
ing a constant bath density of states and g-independent 
coupling Vi : q — V\, the half- width of level 1 due to 
the coupling to the conduction electrons is given by 

Tl = V\ 2 ■ Pbath- 

Now in the situation where the two impurity levels are 
well below the Fermi energy of the conduction electrons 
(cf > 2e + U'), and the intra-level Coulomb repulsion 
U is strong enough to prevent double occupation of each 
impurity level (2e + U > ep) the impurity will be doubly 
occupied. And due to Hund's rule coupling the impurity 
will then be in a total spin-triplet state, i.e. will have 
total spin Si — 1 (See App. [X]for further details). 

In this situation, switching on the coupling Ht be- 
tween the impurity and the conduction electron bath 
gives rise to hopping processes as depicted schematically 
in Fig. [3]which will partially screen the total spin-1 of the 
impurity by flipping the spin in impurity level 1. This 
partial screening of the impurity spin 5 > 1/2 by a sin- 
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gle conduction electron channel is called underscreened 
Kondo effect (UKE)<£2 The coupling to the residual spin 
in the other impurity level gives rise to a so-called sin- 
gular Fermi liquid (SFL) bchaviour 64 i 71 i 72 , in contrast 
to the normal Fermi liquid behaviour of the usual fully 
screened Kondo effect. 

The SFL is characterized by a cusp in the spectral den- 
sity at low temperatures, i.e for low energies the spectral 
density of impurity level 1 is approximately given by : 70 ! 73 



(10) 



where Tq is a new temperature scale associated with the 
Spin-1 UKE and b > is a constant. The cusp in the 
spectral density is related to a logarithmic divergence 
of the quasi-particle weight in the underscreened Kondo 
regime: 65 z oc 1 / u log(kTo / ui) . Hence in contrast to 
the normal Kondo effect there is no well defined quasi- 
particle associated with the UKE. Hence the name sin- 
gular Fermi liquid. Note that although the quasi-particle 
weight z diverges for u 0, the spectral density pi{ui) 
itself does not diverge. 

In an actual experiment the logarithmic cusp charac- 
teristic for the UKE is probably hard to resolve due to 
limited resolution and the smoothening effect of finite 
temperature. Hence in practice the zero-bias anomaly 
arising from the UKE is undistiguishable from that aris- 
ing from the normal Lorentzian-type Kondo peak unless 
very low temperatures can be reached and the experi- 
mental resolution is fine enough to resolve the cusp. 

The half width of the resulting UKE resonance 
is determined by the temperature Tq as: Tf^ 1 = 
kTo exp(— y/2b). In order to compare this width with 
the one of the normal Kondo peak in the 1AIM we 
consider the particle-hole symmetric regime. Following 

Ref. [z3 we have kT oc exp (—£^(U + Jh/2)) while 

kTK oc exp(— (7rf/)/(4Td)). Hence for the same parame- 
ters U and = Tx, the width of the resonance should be 
smaller in the case of the UKE than for the normal Kondo 
effect due to Hund's rule coupling Jh and the additional 
exponential factor exp(— V2b) < 1. Assuming that the 
interaction U and coupling T for the d-level giving rise 
to the Kondo effect is approximately the same for the 
three transition metals considered here, this might ex- 
plain why the widths of the zero-bias anomalies obtained 
for the cases of Fe and Co are considerably smaller than 
for Ni (see Sec. E}. 

Finally, for the 2AIM in the UKE regime one obtains 
exactly the same formula relating the occupation n\ of 
the impurity level 1 (the one coupled to the conduction 
electron bath) to the width and position of the resonance 
as in the case of the 1AIM, eq. ([5]). 70 This explains why 
the results in the case of Co and Fe can be fitted so well 
to the formula for the 1AIM (see Sec. |V|. 



D. Kondo- Fano lineshapes 

The goal of this section is to devise a simple model 
in order to understand the occurrence of different Fano 
lineshaps (i.e. peaks, dips or asymmetric Fano curves) for 
the same material. Our simple model also demonstrates 
the complicated dependence of the Fano parameters on 
the basic microscopic parameters. Fano lineshaps in a 
related model have recently been studied by Zitko 7 ^ us- 
ing the Numerical Renormalization Group and focusing 
on the temperature dependence of the Fano lineshapes. 
Here we neglect any temperature effects since the mea- 
sured Kondo scales are much higher than the tempera- 
ture of 4.2 K at which the experiments where performed. 
In any case the goal of this section is not to give an ex- 
act description of the conductance spectra but rather to 
obtain a qualitative understanding of how different line- 
shapes emerge and how they depend on the microscopic 
details of the system. 

In the following we will show that the low-bias conduc- 
tance (||eU|| < Tk) of our simplified model of a nanocon- 
tact shown in Fig. [5] is well described by the Fano for- 
mula: 



G{V) = g oS + 



A 



1 



q 2 e 2 



(6 + q)2 with e =^^ (11) 



1 



T K 



where g s is the conductance offset, A is the amplitude of 
the Fano resonance and q is the Fano factor determining 
the shape of the Fano resonance (see Fig. |4|), and as 
described before ex is the energy position and Tk the 
half-width of the Kondo peak in the d-level which thus 
determine the position and width of the resulting Fano 
lineshape. 



q ~» 00 / 


dl/dV 


q= 1 

goff ^ 
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FIG. 4. Graphical definition of the Fano resonance parame- 
ters. 

As can be seen from Fig. @] when the Fano factor 
becomes very large (q — > oc) the conductance has a 
Lorentzian lineshape. This is the case when the coupling 
V s d between the s- and the d-level of our tip atom is negli- 
gible so that the zero bias anomaly in the conductance is 
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dominated by the direct transmission through the Kondo 
resonance of the d-level. On the other hand for q = we 
obtain a dip-feature in the conductance. This is for ex- 
ample the case when the direct transmission through the 
d-channel becomes negligible (e.g. for Tr^ ~ 0) so that 
the conductance is only given by the s-channel which 
features a Lorentzian dip due to the coupling V sd to the 
Kondo resonance in the d-level. For |g| = 1 the Fano 
formula gives the typical asymmetric lineshapes. In this 
case the conductance is also dominated by the s-channel 
coupled to the Kondo resonance in the (i-level as for q = 
but now the s-level of the tip atom is not near the Fermi 
level. 

In the appendix we give a derivation of the Fano for- 
mula (llip for our simplified model of a nanocontact 
shown in Fig. [5] The (left) tip atom where the Kondo 
effect is taking place is modeled by one s- and one (i-level. 
While the s-level couples well to both electrodes (via Tl }S 
and r# jS ) and thus has a nearly perfect transmission, the 
d-level hosting the Kondo resonance has a much weaker 
coupling to both electrodes (Tl,<i and T^)). Generally 
we assume that the couplings to the left electrode are 
stronger than to the right electrode, i.e. > T^ a . 

Additionaly there is a small hybridization V s d between 
the s- and the d-level due to the crystal field. 

As Meir and Wingreen showed in their landmark pa- 
per, at zero temperature and in linear response the con- 
ductance through a nanoscopic conductor is well de- 
scribed by the Landauer formula even in the case of a 
strongly interacting system^ The conductance G for 
small bias V is then given in terms of the quantum me- 
chanical transmission function T(u>) as G(V) — Go X 
T(eV) where Go = 2e 2 /h is the fundamental conduc- 
tance quantum. 

As is shown in the appendix, the transmission T(to) 
through the tip atom can be decomposed into the contri- 
butions of direct transmission through the indivdiual s- 
and d-channel, namely T s (uj) and Td(u>), respectively, and 
a mixed channel involving hopping between both chan- 
nels, T sd (uj): 



T(w) = T a (u) + T d (tu) + T sd {w) 



(12) 



Following eqs. (|B8IB9I) the direct channel transmissions 
T a and T d are given by the spectral densities of the s- 
and d-level, respectively, and the couplings of the s- and 
d-levels to both electrodes L and R. 

The Kondo effect gives rise to the appearance of a 
Kondo resonance in the the d-level spectral function 
p d (to) given by eq. ([3]). Hence the contribution of the 
d-channel to the total transmission has a Lorentzian line- 
shape: 



T d {w) 



L,d ■ 1 R,d 



4r, rf -r 



L,d • 1 R,d 



1 



(u-e K ) 2 +T 2 K T LA + T R , d l + x 2 



(13) 



where we have used Tk = z{T^ d + T R d )/2 and we have 
defined the dimensionless quantity x — (to — 6k)/^k- 
Hence as pointed out above, the d-channel contribution 



to the transmission can only give rise to Fano-lineshapes 
with q — > oo. 

We assume that the unperturbed s-channel (i.e. with- 
out coupling to the d-level) has a featureless (i.e. flat) 
and almost perfect transmission T s ° w 1. Due to the 
coupling V sd to the d-level this transmission is modified 
according to (see appendix): 



1 



2zV 2 d ( x + q ) 2 -(x 2 + l) 



T s ■ T K {x 2 + !)•(! + g 2 ) 



(14) 



where T s =Tl,s + ^r,s is the total broadening of the s- 
level and the dimensionless quantity go has been defined 
as the ratio between the s-level energy and the s-level 
broadening: go = — 2e s /r s (see appendix for details). 
The second term of the r.h.s. is of the Fano form, eq. (Tn| 
and represents the modulation of the almost perfectly 
transmitting s-channel due to the coupling to the Kondo 
resonance in the d-channel. Since generally 2|e s | < T s 
we should have |go| < 1 and therefore the s-channel con- 
tribution to the transmission (TT4|) can only give rise to 
dip-like (go ~ 0) or asymmetric Fano lineshapes g ~ 1 
but not to the peak lineshapes where |go| 3> 1. 

For the mixed-channel contribution T sd to the total 
transmission we find the following expression in the ap- 
pendix: 



z 2 V 2 

T sd {u) = T° • f-fa 

1 K 



R,d 
^R,s 



r L,d 

Tl, s 



1 



1 



(15) 



This contribution describes transmission processes where 
an electron hops from one electrode to the s-level of the 
tip atom, subsequently to the d-level via V sd , and then to 
the other electrode. Due to the Kondo peak in the d-level 
it gives rise to a Lorentzian lineshape in the transmission. 

Hence we have shown that our model can give rise to all 
possible Fano lineshapes as obtained in the experiments. 
More specifically, the s-channel contribution T s (w) can 
give rise to the dip-like features (g w 0) and the asym- 
metric Fano features (|g| 1) while the d-channel contri- 
bution T d (uj) and the mixed channel contribution T sd (uj) 
give rise to Lorentzian lineshapes (|g| — ¥ oo) in the trans- 
mission. Which term dominates depends on the specific 
amplitudes of the different transmission channels given 
by the basic parameters of our model. 

It is of course possible to achieve the "canonical" Fano 
lineshape form for the conductance as in eq. pT|) by 
summing up all the individual contributions to the total 
transmission (jT2"j) and reorganizing the terms. We then 
obtain a "new" Fano factor q different from the Fano 
factor go for the pure s-channel contribution (|14[) . This 
new Fano factor will depend on go and the amplitudes 
of the individual contributions to the transmission, and 
yields a relatively complicated expression in terms of the 
basic parameters of our model. The same is true for the 
amplitude A of the Fano feature defined by eq. (fTTj) . 

However, we can obtain quite simple expressions for 
g and A in eq. (1111) in an important limit of our model, 
namely when the coupling of the d-level giving rise to the 
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Kondo peak to one of the electrodes becomes very small, 
e.g. Tjid — > 0. In that case the direct transmission 
through the d-channel is strongly suppressed, i.e. TjaO 
so that now only the mixed channel contribution T a d can 
give rise to a Lorentzian lineshape in the transmission. 

q = \ F" ■ q ° = — /F 2£ p (16) 
i + r a /r L , a 4V% 

1 + ql Y s ■ T L4 

'^■k. T ^- v -'- AK ,17 » 

where Ak is the amplitude of the Kondo resonance in the 
spectral function of the d- level as given by eq. (j8|) . Note 
that the Fano factor qo of the s-channel transmission T s 
is now scaled by the factor y/T s /TL^ s > 1 to yield the 
Fano factor q of the Fano lineshape in eq. ([TTj) meaning 
that for large values of the ratio -\/T s /Tl -s we can obtain 
\q\ values > 1 for < |go| < 1- Furthermore, we see that 
the amplitude of the Fano resoance in the conductance 
is proportional to the amplitude of the Kondo resonance 
and to the square of the coupling V s d between the s-level 
and the <i-level of the tip atom. 

IV. EXPERIMENTAL DETAILS 

The experiments were performed using a home-made 
Scanning Tunnelling Microscope (STM) operated in a He 
cryostat at 4.2K. Two pieces of the same metal wire (Fe, 
Co or Ni) of 0.1 mm of diameter were scratched and son- 
icated in acetone and isopropanol before being mounted 
as 'tip and sample' in the microscope. The conductance 
between the two pieces of metal is obtained in a two ter- 
minal configuration by measuring the current at a fixed 
bias voltage, in this case 100 mV. In these conditions we 
can record traces of conductance wile changing the dis- 
tance between the two metals (a typical trace is shown in 
Fig. Hk) in a similar fashion as performed in other break 
junction experiments^. The samples are then prepared 
at low temperatures by indentation until no traces with 
subquantum events are shown. The histograms (Fig. [T]d) 
are similar to those at Mechanically Controlled Break 
Junction (MCBJ) experiments where a fresh surface is 
formed at cryogenic vacuum when breaking a notched 
wire by the controlled bending of a elastic substrate (See 
e.g. Ref. for details). This shows that our mea- 
surements are performed over a clean spot of our sam- 
ples. A strong indentation between the two electrodes is 
performed between the fabrication of consecutive atomic 
contacts to ensure the cleanliness. Our STM set-up has 
imaging capabilities, however, our surfaces are not atom- 
ically flat due to preparation and in order to acquire a 
large number of contacts to analyse no imaging is per- 
formed between contacts formation. 

As described in Sec HH traces show plateaus coming 
from the atomic rearrangement of the wire while pulling 



and the last plateau around the quantum of conductance 
is associated with the formation of a single atom contact. 
Once the conductance of the one-atom contact has been 
determined by the position of the first peak of the con- 
ductance histogram, we can fabricate monoatomic con- 
tacts by stopping the breaking process of the contact at 
the desired value of conductance within the range de- 
fined by the conductance histogram. The stability of 
our system allows us to maintain such an atomic con- 
tact for hours. We study the transport spectroscopy of 
these contacts in a similar way as performed in other 
transport experiments either in tunnelling or high con- 
ductance regimes (e.g. the case of quantum dot devices 
or the spectroscopy of adatoms by STM). We sweep the 
bias voltage from -100 to 100 mV while recording the 
dl/dV signal with the help of a lock-in amplifier when 
adding a lmV AC excitation at a frequency about lKHz 
to the applied bias voltage. 



V. RESULTS 

The fabrication of atomic contacts by using a STM (or 
MCBJ) offers the possibility of studying a high number of 
different contact configurations in a reasonable amount 
of time. We fabricated hundreds of atomic contacts of 
Fe, Co and Ni and performed electron spectroscopy mea- 
surements as described in Section llVl About 80% of these 
curves showed clear asymmetric profiles centred at zero 
bias as the ones plotted in Fig. [SJ 

Similar asymmetric zero bias anomalies (ZBA) were re- 
ported for ferromagnetic atomic contacts and attributed 
to the existence of conductance fluctuations 76 or to the 
existence of a magnetic domain wall. 77 On the other 
hand, these asymmetric profiles resemble the data re- 
ported for single magnetic adatoms in the contact regime, 
not only in shape but also in the energy scale of the 
features^. 

In our previous work we have shown clear evidence 
of the Kondo effect being responsible for these ZBA. 38 
Recently, ZBA profiles, possibly of similar origin, have 
been reported in Refs. [78j and l79l for nanocontacts made 
from other materials and also show the Kondo related 
ZBA in adatoms contacted by ferromagnetic tip o 36 ' 37 . 

As in other Kondo systems, the asymmetric lineshapes 
can be fitted to the Fano equation (| 1 1 [) . From this fitting 
we extract the values for the different parameters that 
describe the Kondo effect. As sketched in Fig. [H the 
width of the resonance is directly related to the Kondo 
temperature (2V). The Fano parameter q contains infor- 
mation about the symmetry of the lineshape. We denote 
by ex the energy at which the resonance is centered. As 
introduced in Sec. IIII1 this parameter is associated with 
the energy position of the Kondo resonance and therefore 
to the occupation of impurity level rid, eq. ([6]). Finally, 
A is the amplitude of the Fano profile and g g the con- 
ductance offset of the curve out of the resonance. 

Each realization of the contact leads to a slightly differ- 
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FIG. 5. Representative differential conductance curves for 
atomic contacts of Fe, Co and Ni. A characteristic resonance 
associated with the Kondo effect appears at small bias. The 
corresponding fit to the Fano lineshape is shown in black for 
each of these resonances. Due to the large differences in the 
width of the resonances, the fitting range is chosen so that is 
possible to fit the whole range of the resonance to the Fano 
lineshape. The selected curves exemplify different symmetry 
cases, from the more symmetric {q ~ oo and q ~ respec- 
tively) to the clearly asymmetric ones (q ~ ±1). While all 
possible symmetries seem to happen for the three materials, 
the width of the resonance which is associated with the Kondo 
temperature, is different for the three materials, especially for 
the case of Ni. 

ent configuration. The statistical distributions of the dif- 
ferent parameters of the Fano equation described above 
will reflect the subtle differences in the electronic struc- 
ture of each contact. We present below a statistical anal- 
ysis of each of these parameters for hundreds of contacts 
of Co, Fe and Ni. This novel statistical analysis (since in 
our previous work2£ we only analysed briefly the shape 
and mean value of the Kondo Temperature) together 
with theoretical considerations brings new insight into 
the physics of Fano-Kondo resonances. We present now 
the distributions of these paramctercs and compare the 
results between materials. 



A. Kondo Temperatures 

As described in Sec. IIIII the zero-bias resonances ob- 
served in the conductance characteristics of Kondo sys- 
tems are directly related to the resonances developed in 
the spectral density of the system. As explained in the 
theory section, the width of this resonance is determined 
by the energy scale of the Kondo screening, the so-called 
Kondo temperature Tk ■ Thus the width of the observed 
Fano lineshapes must be proportional to the Kondo tem- 
perature of the system. More precisely, we define the 
Kondo scale as the half width of the Fano lineshape: 
Tif = ^bTk- The width of the measured Kondo reso- 



nance is strongly affected by a finite temperature of the 
system: in addition to the standard thermal broadening 
of any differential conductance feature, the Kondo res- 
onance presents an intrinsic thermal broadening^ This 
results in a considerable extra broadening of the reso- 
nance at temperatures on the order of magnitude of Tk- 
In our case, since the width of our resonances excesses 
in more than an order of magnitude the experimental 
temperature of 4.2K and the bias voltages used are low 
enough 3 - , we can disregard thermal effects and consider 
that we can extract the Kondo temperature for each con- 
tact directly from the width of the Fano resonance. 

As we have already described in Ref. HH, the distri- 
bution of Kondo temperatures fits a logarithmic normal 
distribution for the three materials Co, Fe and Ni, mean- 
ing that the logarithm of Tk is normally distributed 
(presented in Fig. [6]). This peculiar behavior is easily 
understood when interpreted in terms of the Kondo ef- 
fect: Since many different atomic configurations result 
in single-atom contacts, their electronic properties, such 
as conductance (Fig. [T]), density of states and the as- 
sociated energy scales are expected to be normally dis- 
tributed. On the other hand following eq. ^ the Kondo 
temperature depends exponentially on the typical energy 
scales of the problem. Hence In Tk for different contacts 
should follow a normal distribution if the relevant energy 
scales of the problem are normally distributed. 

By just looking at the resonances, as for example the 
ones shown in Fig. [SJ it can be observed that the Fano 
features in the conductance spectra of Ni contacts are 
considerably broader than the ones for the case of Co and 
Fe. As shown in Fig. [6] and summarized in table HI the 
histograms yield most frequent values for the resonance 
widths of T K = 90 K, 120 K and 280 K for Fe, Co and Ni 
respectively, following the same trend T^ e < T^° < T^ 1 
as in the case of adatoms of these elements deposited 
on non-magnetic surfaces^ and diluted alloys of Cu con- 
taining the same concentration of magnetic atoms^. In 
simple terms, the Kondo temperature decreases as the 
size of the screened magnetic moment increases, as we go 
from Ni to Fe. 





Fe 


Co 


Ni 


Tk (K) in bulk Cu^ 


10-50 


300-700 


~1000 


Tk (K) adatoms 




53-9222 


120^ 


T K (K) this work 


90 


120 


280 


m a tom(/i-B) 


3 


2 


1 


m hu]ii (uB) 


2.22 


1.72 


0.60 



TABLE I. A comparison of Kondo temperatures in the case 
of magnetic impurities in bulk Cu, magnetic adatoms in tun- 
nelling regime and for the ferromagnetic contacts (Fe, Co and 
Ni) in this work. The different values of Kondo temperatures 
in the case of adatoms correspond to measurements performed 
over different substrates. Also shown are the magnetic mo- 
ments of isolated atoms and in bulk. 

Figure |5] shows the distribution for ln(Tff) fitted to a 
Gaussian distribution. Surprisingly enough, in spite of 
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Ln(T K ) 
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Fe 


Co 


Ni 


Mean Tk 


90 


120 


280 


FWHM In(Tk) 


0.96 


0.91 


0.73 
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FIG. 6. Histograms of values of logarithm of Kondo temper- 
ature for Fe, Co and Ni. The histograms are normalized to 
the total number of counts. The distribution of values for the 
ln(Tk) is clearly narrower for Ni than for Fe and Co in spite 
of the higher value of Kondo temperatures for Ni, as sum- 
marized in the inset. The most frequent value of the Kondo 
temperature is determined with an error of 10 K. 



showing quite higher values of Kondo temperatures, the 
distribution of Ni when plotted in logarithmic scale is 
considerably narrower than in the case of Fe and Co (see 
also inset of Fig. This suggests that in the case of 
Ni the characteristics of the Kondo screening are very 
different from the cases of Co and Fe possibly indicating 
a different mechanism for the case of Ni and the cases of 
Co and Fe. 

The higher Kondo temperatures for Ni as well as their 
narrower distribution could well be connected to the dif- 
ferent chemical valence and the resulting magnetic mo- 
ment in comparison to Co and Fe: While Ni basically 
has one hole in the 3d-shell and therefore features an 
atomic spin of 1/2, Co and Fe have two and three holes 
in their 3d-shells associated with atomic spins of 1 and 
3/2, respectively. Hence in the case of Co and Fe the pos- 
sibility exists that the full atomic spin S > 1/2 is only 
partially screened while the spin-1/2 in the case of Ni is 
fully screened by the conduction electrons. The cases of 
Co and Fe would then resemble the situation of an under- 
screened Kondo lattice where the remaining unscreened 
spin couples ferromagnetically to the spins on neighbour- 
ing atoms as discussed in Ref. l63l As explained in the 
theory section IIIII such an underscreened Kondo effect is 
characterized by sharper resonances in comparison to the 
normal fully screened Kondo situation, and hence results 
in lower Kondo temperatures in the analysis. In this 
sense the Kondo temperatures may be underestimated 
in the case of Co and Fe. We discuss this in detail in 
the discussion section IVTl together with the distribution 
of other parameters. 



B. Resonance energy (ek) and d- level occupation 

Another important parameter to study is the position 
of the Kondo resonance €k which accounts for the devi- 
ation of the center of the Fano resonance from zero bias 
and is related to the energy of the effective Kondo level 
with respect to the Fermi energy^. The inset of Fig. [7J 
shows the distribution of this parameter for hundred of 
contacts of the three materials under study. These dis- 
tributions fit a Gaussian lineshape which is considerably 
broader for the case of Ni than for Co and Fe. 
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d-level Occupation (n d ) 



FIG. 7. Histograms of the d-level occupations nd for contacts 
of Fe, Co and Ni. The color lines show the fitting of these dis- 
tributions to a Gaussian peak. Surprisingly, the distribution 
for Ni is narrower than the other two, in spite of the fact of 
this material showing a broader distribution of the parameter 
en- Inset: Distribution of en for hundred of contacts of Fe, 
Co and Ni. This parameter accounts for the deviation of the 
center of the Fano resonance from zero bias and is related to 
the position of the localized magnetic moment in the Kondo 
model. 

Most interestingly, from the values of €k and Tk it is 
possible to extract the occupation of the <i-level (giving 
rise to the Kondo resonance) from the experimental data, 
assuming a Fermi liquid approximation as explained in 
Sec lIH( eq. (|6|). This approximation should be valid here 
since we are well below the Kondo temperature of our 
system. As explained in Sec. IIIII this relation is even 
valid for the underscreened Kondo effect in the multi- 
orbital Anderson model (although strictly speaking we 
then have a Singular Fermi liquid) where now rid refers to 
the the d-level whose spin is screened by the conduction 
electrons, and ek and Tk refer to the position and width 
of the Kondo resonance of the underscreened Kondo ef- 
fect. Fig. [7] shows the distribution of positions and the 
resulting calculated occupations nd for Fe, Co and Ni 
which again follow a normal distribution. 

As summarized in Table Hi] the mean value for the 
three materials is close to 1. Interestingly, in spite of 



the broader distribution of values of €k, the occupations 
for Ni contacts clearly show a narrower distribution. Nu- 
merically this can be explained by the much narrower 
distribution of Itl(Tk)- Physically, the reason behind the 
narrower distribution of Kondo temperatures and occu- 
pations in the case of Ni might be that charge fluctuations 
become stronger with an increasing number of active lev- 
els as already discussed in the seminal work of Nozieres: 62 
In the case of Ni we are most likely dealing with a single 
active impurity level due to a single hole in the 3d-shell 
of Ni. Hence the charge in this level is quite well defined, 
and the occupation very close to one electron. In the 
cases of Co and Fe on the other hand we should have 
more than one active impurity level and therefore the 
variation in the occupation of the level giving rise to the 
Kondo resonance is much stronger. Relatedly, a broader 
distribution in the occupation of the impurity level giv- 
ing rise to the Kondo resonance should also give rise to 
a broader distribution in the Kondo temperatures since 
charge fluctuations strongly alter the width of the Kondo 
resonance £ 





Fe 


Co 


Ni 


Mean ex (mV) 


-1.0 


0.6 


0.3 


FWHM ex(mV) 


14.7 


16.3 


25.4 


Mean rid 


1.07 


0.97 


0.99 


FWHM n d 


0.79 


0.74 


0.52 



TABLE II. Most frequent value of Resonance energy and d- 
level occupations and their respective width of their Gaussian 
distributions for Fe, Co and Ni 



C. Amplitudes 

At first glance, the distribution of amplitudes in Fig. 
[8] shows no clear differences between the three materi- 
als. The most frequent value of the amplitude is approx- 
imately 0.1 2e 2 /h, i.e. about 10 percent of the conduc- 
tance of the contact. On the other hand we find that 
the square root of the amplitude indeed follows a nor- 
mal distribution as can be seen from the inset of Fig. 
[5] which shows the statistical distributions of the square 
root of the amplitude for the three materials. A possible 
explanation is that following expression (|17[) the ampli- 
tude depends quadratically on the coupling V s d between 
the s- and cZ-level of the tip atom. V s d is expected to 
vary strongly when the atomic configuration of the con- 
tact changes since it is induced by disorder in the contact 
region, and is absent for perfect crystalline order. Hence 
if V sc i is normally distributed and is the parameter de- 
termining A that is most strongly affected by changing 
the atomic configuration of the contact one would expect 
\[A~ to be normally distributed. 

Furthermore one can seen that the distributions of \J~A~ 
for the three materials are centred at quite similar values 
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Amplitude (2e 2 /h) 

FIG. 8. Distribution of amplitudes extracted from the fitting 
of characteristics of hundred of contacts of Fe, Co and Ni to 
the Fano equation. The distributions are similar for the three 
materials, being the amplitudes about a 10 percent of the con- 
ductance of the contacts. The inset shows the distribution of 
and its fit to a Gaussian. This distribution could reflect 
the quadratic dependence of amplitude in different coupling 
terms. 

and also have similar widths. However, we can make 
out a subtle trend (Tab. Illlf) : The mean amplitude and 
the mean value of the square root of the amplitude both 
are slightly higher for Fe contacts than for the other two 
materials, and also both distributions are slightly broader 
for Fe than for Co and Ni. However, this trend is not as 
clear as the trend observed in the distributions of Kondo 
temperatures Tk for the three materials (Fig. [6]). Thus 
it is difficult to draw any further conclusions from it. 
Moreover, we would like to point out that a similar trend 
is observed in the distributions of the conductances (see 
Fig. [T]and Tab. |Hl|: Fe has a higher average conductance 
than Ni and Co. More data would be needed to extract 
further conclusions from this analysis. 





Fe 


Co 


Ni 


Mean A from \[A 


0.13 


0.088 


0.10 


Mean VA 


0.36 


0.30 


0.32 


FWHM VA 


0.31 


0.24 


0.27 


Mean G from Fig. [[] 


2.0 


1.2 


1.6 



TABLE III. Mean values of amplitudes and width of the the 
gaussian distributions of \J~A and mean conductances for Fe, 
Co and Ni. 

Far more interestingly, the plot of the amplitude versus 
the Kondo temperature presents an intriguing trend as 
can be seen from Fig. [9] For all three materials there 
seems to be an approximately linear dependence of the 
amplitude on the Kondo temperature with respect to the 
Kondo temperature. Moreover, when normalizing the 
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amplitudes and Kondo temperatures to the respective 
mean values (Tab lIII[ ) one obtains a universal dependence 
suggesting that A oc f(Tx)- 

This universal behaviour seems to reflect some kind 
of universality in the relation between the basic param- 
eters U, ed and Td that ultimately determine the Kondo 
properties of our system. Said in another way, the basic 
parameters U , td and are not independent from each 
other but are linked together in such a way that univer- 
sal scaling between the amplitude Ak and the Kondo 
temperature Tk results. For example it is conceivable 
that both Yd (the coupling of the impurity level to the 
conduction electrons) and U (the effective Coulomb re- 
pulsion of the impurity level) are related since a change 
in the coupling Td implies a change in the localization of 
electrons in the impurity level and hence can result in an 
alteration of the screening of the effective Coulomb inter- 
action U. Further theoretical work is necessary in order 
to achieve a rigorous interpretation of these results. 
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FIG. 9. Left: Scattered plot of the amplitude of the reso- 
nance versus the Kondo temperature obtained for hundreds 
of contacts of Fe, Co and Ni. The plot shows a clearly similar 
(linear) trend for the three materials. When divided by the 
most probable value of the Kondo temperature for each ma- 
terial, the distributions lay over each other, as shown in the 
right panel. 



D. Fano Parameters 

The Fano parameter q accounts for the symmetry of 
the Fano resonances, recovering the perfect Lorentzian 
shape for q — > oo and its inverse for q = 0. As already 
commented above all possible symmetries are found for 
each material. A histogram of the values of q for the 
three materials is shown in Fig. I10I Since q ranges from 
to oo, an alternative representation where q — tan(a) 
is chosen for simplicity. In this representation a = cor- 
responds to the dip and a = ir to the peak lineshapes. 
The totally asymmetric cases are represented by a = n/2 



(q = 1) and a = 3tt/2 (q = —1). There is a remarkable 
preference towards the symmetric cases over the more 
asymmetric ones in the case of Fe and Co, and in par- 
ticular the dip-like (q = 0) lineshapes seem to prevail. 
For Ni, the occurrence of asymmetric cases (\q\ ~ 1) is 
considerably higher, overall the Fano factors seem to be 
more evenly distributed than in the case of Fe and Co. 
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FIG. 10. The so-called Fano parameter q ranges from to oo, 
but it can be expressed as q = tan(a) with a between and n, 
for simplicity in the representation, a = corresponds to the 
deep, a — n to the peak and a = 1, —1 the totally asymmetric 
cases. The histograms are normalized to the total number of 
data. 

In the case of adatoms in tunnelling, the symmetry is 
well understood in terms of the ratio between the proba- 
bility of transmission through the d or the s channels^. 
Nevertheless, in the point contact regime and in our case, 
the system presents a higher complexity. There is a to- 
tally open s-channel and the nature of the orbitals and 
couplings contributing to the interference is not com- 
pletely clear. The distribution of the q parameters should 
contain information about the contributions from the dif- 
ferent channels to the interference which might be ex- 
tracted with the help of an appropriate theoretical treat- 
ment. 

In Fig. [TT] we plot the evolution of a single Co contact 
when stretching the junction. The symmetry of the curve 
changes as the contact rearranges in a different manner 
which translates into a different conductance of the curve 
(the curves are not offset). When stretching the con- 
tact the the hoppings between different orbitals, the elec- 
tronic structure and hence the conductance varies, this 
affects the interference between the conduction channels 
and also the Kondo resonance responsible for the Fano 
lineshapes. This evolution is probably specific for each 
contact since an attempt to find a statistical relation be- 
tween q and the conductance offset results in a random 
distribution of points. More work is needed to fully clar- 
ify this effect. 
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FIG. 11. Evolution of the dl/dV characteristics of a single Co 
contact as stretching the junction. As the conductance of the 
contact decreases, also the symmetry of the curve, that is the 
Fano parameter, q changes as a consequence of the different 
electronic structure of different arrangements. 



E. Summary 

In summary, the statistical analysis of the Fano-Kondo 
lineshapes extracted from hundreds of atomic contacts 
shows that: 

(i) The statistical distribution of Kondo temperatures 
for each material is a log-normal distribution, i. e. the 
logarithm of the Kondo temperatures is normally dis- 
tributed. This is expected in the context of the Kondo ef- 
fect since the Kondo temperature depends exponentially 
on the different electronic properties and couplings of the 
system. These properties are likely to be normally dis- 
tributed for the slightly different atomic configurations 
of the contacts. Hence the log-normal distribution of 
Kondo temperature presents further evidence that the 
Kondo effect is responsible for the zero-bias anomalies in 
the conductance. 

(ii) The mean value for the Kondo temperatures is con- 
siderably higher for Ni than for Co and Fe whereas the 
distribution of the logarithm of the Kondo temperature 
for Ni is narrower than for Co or Fe. This finding sug- 
gests that possibly a different mechanism is responsible 
for the Kondo screening in the case of Ni than in the case 
of Co and Fe. 

(iii) The distribution of the resonance energies ex is 
Gaussian and is centred around zero for all three mate- 
rials. In the case of Ni this distribution is much broader 
than in the case of Fe and Co where the distribution 
widths are similar. This again points to different mech- 
anisms for the Kondo screening for Ni on the one hand 
and for Co and Fe on the other hand. 

(iv) The occupations of the impurity level giving 
rise to the Kondo resonance calculated from Tk and ck 
by eq. ^ again follow a Gaussian distribution centered 



around occupation 1. This distribution is narrower for Ni 
than for Co and Fe. This can be understood by consid- 
ering that most likely for Ni only a single d-level is active 
while for Fe and Co several d-levels must be active. This 
leads to stronger charge fluctuations and hence a larger 
variation in the d-level occupations in the case of Co and 
Fe. 

(v) The distribution of the Fano curve amplitudes does 
not show significant differences between the materials. 
We find that the square root of the amplitudes is nor- 
mally distributed. This again can be understood by 
a normal distribution of the characteristic parameters 
of the nanocontacts on which the amplitude depends 
quadratically. 

(vi) The plot of the amplitudes against the Kondo tem- 
peratures follows a similar (almost linear) trend for the 
three materials. When divided by the mean values for 
each material the data lines are perfectly on top of each 
other, showing a universal scaling behaviour of the am- 
plitude with the Kondo temperature. 

(vii) There are more asymmetric Fano lineshapes in the 
case of Ni than for Co and Fe. The distribution of Fano 
parameters q again is similar for Fe and Co but different 
for Ni. Co and Fe show a clear preference for the dip-like 
q = line shapes while for Ni the distribution of the Fano 
parameters is somewhat more uniform. 

Table IIVI shows a summary of the parameters ex- 
tracted from the statistical analysis of the Fano-Kondo 
lineshapes for the all three materials. 





Fe 


Co 


Ni 


Mean T K 


90 


120 


280 


FWHM ln{T K ) 


0.96 


0.91 


0.73 


Mean ck (mV) 


-1.0 


0.6 


0.3 


FWHM ex(mV) 


14.7 


16.3 


25.4 


Mean rid 


1.07 


0.97 


0.99 


FWHM n d 


0.79 


0.74 


0.52 


Mean A from </A 


0.13 


0.088 


0.1 


FWHM yT(A) 


0.31 


0.24 


0.27 



TABLE IV. Summary of the different parameters extracted 
from the statistical analysis of hundreds of atomic contacts 
for each of the three materials Fe, Co and Ni. 



VI. DISCUSSION 

The picture of the Kondo effect in ferromagnetic 
atomic contacts that emerges from our statistical anal- 
ysis and theoretical considerations is the following: The 
low coordination and disorder of the atoms in the con- 
tact region in connection with a higher effective Coulomb 
repulsion can lead to the localization of a single spin in 
an individual cZ-level of an atom in the contact region.— 
Now, due to disorder in the contact region it can happen 
that this ci-level only couples very weakly to the spin- 
polarized e?-levels on neighbouring atoms but instead it 
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has an effective coupling to the basically spin-unpolarized 
s-type conduction channel. In this situation the Kondo 
effect can take place and screen the spin in that d-level. 

In the case of Ni there is one hole in the 3d-shell and 
therefore a spin- 1/2 associated with it. If the spin be- 
comes localized in a <i-level that predominantly couples 
to the spin-unpolarized s-type conduction electrons this 
spin can be fully screened. Hence for Ni we should have a 
normal Kondo effect where the full spin-1/2 is screened. 
In the case of Fe and Co the atomic spin is higher due to 
3 and 2 holes, respectively, in the 3d-shell of these atoms. 
This might explain why the measured Kondo tempera- 
tures are higher in the case of Ni than for Co and Fe since 
the Kondo temperature decreases with increasing spin of 
the impurityj 81 i 82 However, one would then also expect 
a significantly lower Kondo temperature for Fe than for 
Co which is not the case. 

The scenario that we propose instead for Co and Fe 
is an underscreened Kondo effect where the full atomic 
spin S > 1/2 is only partially screened by the conduction 
electrons. More precisely, only the spin-1/2 within the 
d-level that predominantly couples to the s-type conduc- 
tion electrons will be screened while the rest of the spin 
S — 1/2 which is likely to be localized in d- levels that 
couple more strongly to the spin-polarized d-levels on 
neighbouring atoms, remains unscreened. This scenario 
fits very well with our results: In particular, it explains 
why the average Kondo temperature is very similar for 
Fe and Co but is considerably higher for Ni as the un- 
derscreened Kondo effect is characterized by a cusp-like 
resonance much sharper than the Lorentzian-type reso- 
nance of the normal fully screened Kondo effect. An un- 
derscreened Kondo effect for Fe and Co contacts would 
also explain why the occupations can be calculated from 
eq. (J6j> derived for the single-level Anderson impurity 
model. 

A recent work by Neel et a£l shows how the Kondo 
temperature increases for a Cobalt impurity when it is 
contacted by an Fe tip compared to the case of using a 
Cu tip. These results and their interpretation are in good 
agreement with the results presented here: differences in 
the composition and geometries of the junctions lead to 
different electronic structure of the contact which deter- 
mine the Kondo screening. By studying a specific system 
with a well characterized geometry, Neel et al. can ex- 
tract the change in hybridization and the corresponding 
occupation of the d-level as the main cause of the ob- 
served increase in Kondo temperatures. We study here 
instead a large variety of possibilities. Changes in the 
hybridization are probably responsible for the broad dis- 
tribution of Kondo temperatures in our data, but other 
parameters such as the occupancy of the 3d-shell may 
also change in our system from contact to contact. 

Finally, we would like to discuss the connection to a 
related work by Bork et al^ with the results presented 
here: In the work of Bork et al. the Fano-Kondo line- 
shapes of a system consisting of two Co atoms, one on 
a Cu surface and the other attached to Cu STM tip are 



recorded while the distance between the two Co atoms 
is decreased from the tunnelling regime to the contact 
regime. The authors report a splitting of the Fano-Kondo 
lineshapes when the contact regime is entered due to 
the formation of a spin-singlet state between the two Co 
atoms. 

In some cases the Fano-Kondo lineshapes of the fer- 
romagnetic contacts measured here could also be inter- 
preted as a splitting of the Kondo resonance. For exam- 
ple, the evolution of the Fano lineshape of a Co nanocon- 
tact being stretched (Fig. [TT]) possibly shows a small 
splitting for the curve with the lowest conductance. How- 
ever, if this specific feature really shows a splitting of the 
Fano-Kondo resonance it actually occurs in the opposite 
direction as for the system studied by Bork et al., i.e. 
it occurs when the nanocontact is pulled apart and not 
when the tip atoms are brought together as in the case of 
Bork et al. In any case we have not observed such curves 
that could be interpreted as splittings very frequently, 
and the focus of this work is a statistical analysis of a 
large number of different configurations of atomic con- 
tacts while the work of Bork et al. focuses on a very 
specific system. On the other hand we also would like 
to point out that conductance oscillations in the spec- 
troscopy at the atomic scaled further complicate the in- 
terpretation of this kind of curves. 

Still one might wonder why splitting of the Fano- 
Kondo lineshapes in the atomic contacts studied here are 
not observed more frequently. The reason might be that 
indeed the d-level on a tip atom giving rise to the Kondo 
effect just does not couple very well to the d-levels of the 
other tip atom due to the disorder in the contact region. 
Hence the splitting is probably very weak compared to 
the Kondo temperature and therefore not observed, as 
also shown in Refill, where the Kondo resonance remains 
unsplitted even when the Co adatom is contacted by a 
ferromagnetic tip. Another possible explanation could 
be that the weak coupling between the (i-levels giving 
rise to the Kondo effect localized on different tip atoms 
is compensated by the weak spin polarization of the con- 
duction electrons similar to the case of a double quantum 
dot coupled to ferromagnetic leads ^ 



VII. CONCLUSIONS 

Atomic contacts are a unique system to understand 
the dramatic consequences of low coordination and dis- 
order on the electronic transport and magnetic properties 
at the nanoscale. The emergence of the Kondo effect in 
atomic contacts is a good example for this<2& Here we 
have extended our previous work in several aspects in 
order to gain new insights into the nature of the Kondo 
effect in ferromagnetic one-atom contacts: We have pre- 
sented an exhaustive statistical analysis of the Fano res- 
onances in the spectroscopy of ferromagnetic one-atom 
contacts. In particular, we have analysed not only the 
distribution of Kondo temperatures, but also of the res- 
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onance energies, amplitudes and Fano parameters. Such 
an exhaustive statistical analysis has never been per- 
formed before for any Kondo system, and it provides us 
with new information on the nature of the Kondo effect 
in these system. 

From this analysis we have for example obtained the 
distribution of d-level occupations, and the dependence 
of the amplitudes of the Fano resonance on the Kondo 
temperature. We also point out the widths of the pa- 
rameter distributions as an indicator on the robustness 
of the Kondo screening in one-atom contacts of the dif- 
ferent materials. Additionally, we have further developed 
the theory in order to explain the observed qualitative 
differences between Ni on the one hand and Co and Fe 
on the other hand. Furthermore, we have devised a sim- 
ple microscopic model which allows us to understand the 
occurrence of different Fano lineshapes for the same ma- 
terial in terms of variations of the microscopic interac- 
tions. 

Our statistical analysis shows clear differences in the 
Kondo characteristics of different ferromagnetic materi- 
als, i.e. of Fe and Co versus Ni. These differences can 
be explained by the different valences, leading to dif- 
ferent Kondo scenarios. Although this might not be a 
unique explanation for the observed phenomenology, the 
reported results seem to fit well with an underscreened 
Kondo effect for the case of Fe and Co, and a standard 
spin-1/2 Kondo model for the case of Nickel. 

While the different screening conditions might account 
for the differences in Kondo temperatures and occupa- 
tion, the complexity of atomic contacts, especially of 
transitions metals, where many different channels con- 
tribute to the transport, makes it difficult to fully un- 
derstand the details for each single contact. A sim- 
ple interpretation for certain parameters as the exact 
shape of the curves (Fano parameter) or the amplitude of 
the resonances is difficult to achieve as the spectroscopy 
might contain also information about other phenomena 
such as quantum oscillations due to interferences between 
non resonant conduction channels, inelastic processes as 
phonon excitations etc. 

In the recent past, the combination of molecular 
dynamics simulations and ab initio transport calcula- 
tions and experiment has proven to be a successful 
path for understanding the atomic, electronic and mag- 
netic stucture as well as the transport properties of 
nanocontactS i 39 i 50 ' 59 Further work in this direction in 
combination with more sophisticated many-body tech- 
niques capable of describing the Kondo effect such as the 
Dynamical Mean-Field Theor y 60 i 85 should contribute to 
fully understand Kondo physics and the presence of mag- 
netism in these systems and in atomic-size structures in 
general. 
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Appendix A: Two-level Anderson impurity 
eigenstates 

Here we briefly discuss the many-body eigenstates of 
the isolated two-level Anderson impurity model and the 
fluctuations between them that lead to the underscreened 
Kondo effect as discussed in Sec. IIII CI 

The isolated impurity Hamiltonian H\ of the 2AIM 
© is easily diagonalized. The empty impurity state is 
denoted by and has energy Eq = 0. For single oc- 
cupation the four eigenstates and energies are trivially 
given by 

\d 1 ;ia)=dl\0),E 1 = (ia\H T \ia) =e (Al) 

For double occupation, the eigenstates can be separated 
into spin-singlet (total impurity spin S = 0) and spin- 
triplet (5=1) states. The latter are given by 

f|t,t) (M = l) 
\d 2 ;S=l;M) = \l,2) ®l |t,|) + (M = 0) (A2) 

(M = -l) 

where M is the projection of the total impurity spin onto 
the spin quantization axis. The triplet states are all de- 
generate with energy Et = 2e + U' — Jh/4. The spin 
singlet states can be written as: 

\d 2 ;S = 0;i 3 ) = \i,j) + (A3) 

The corresponding eigen energies are E s 2 = 2e+U'+j Jh 
and El, 1 = Ef = 2e + U. Typically, U U' + J H , and 
therefore Eg 1 = Eg 2 > Eg 2 . Due to Hund's rule coupling 
the triplet states are lower in energy by an amount E s 2 — 
E)t = Jh- In Fig. [12] we show the corresponding energy 
level diagram for the 2AIM. U and U' are assumed to be 
big enough to prevent triple and full occupation of the 
impurity. Hence the triplet states comprise the ground 
state manifold of the isolated impurity system. 

Now in the situation of the underscreened Kondo ef- 
fect where only one of the impurity levels is coupled to 
the conduction electron bath, hopping processes between 
the impurity and the bath can only lead to a spin flip 
in the impurity level that is coupled to the bath as illus- 
trated in Fig. |31 These spin flip processes can only lead 
to fluctuations between states with M = 1 and M = 
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FIG. 12. Schematic energy level diagram for the (isolated) 
2-level impurity described by Hi in eq. (|9jl. 



and between states with M = —1 and M = 0. For 



example: \S = 1; M = 1) 



|lt,2t) ( 



|2f) 



|U,2t> H ^ |2t) |lt,2t> (_) - Hence the to- 

tal spin of the impurity is only partially screened by 
these processes. This is the essence of the underscreened 
Kondo effect. 



Appendix B: Derivation of the Fano formula for 
model contact 



We assume the simplified model of a nanocontact de- 
scribed in Sec. IIIII and shown in Fig. [5J We concentrate 
on one of the tip atoms of the nanocontact which in our 
model consists of an (almost) perfectly transmitting s- 
level and the d-level where the Kondo effect takes place. 
Hence the Greens function (GF) of the tip atom can be 
written as: 



G s G s d 

G s d Gd 



(Bl) 



Gd is the Green's function of the d-level which in the 
Kondo regime is given by eq. JH) and thus yields a Kondo 
resonance in the corresponding spectral function as de- 
scribed by eq. ^ at energy €k, half width and the 
quasi particle weight z. 

The unperturbed s-level GF (i.e. without coupling V s d 
to the d- level) is given by 



1 



1 



•»T./2 



ir s /2 



(B2) 



where e s is the energy of the s-level and T s — T^.s + 
Tji^s the width due to coupling to both electrodes. Since 
the energies we are interested in are on the order of the 
Kondo scale, \uj\ ps Tk and Tr *C T s , we neglect the 
w-dependence of G® in the last step of eq. (|B2I) . 

Now due to the coupling V s d of the s-level to the d- 
level, the full GF G s of the s-level is modified according 
to: 

G s (oj) = G2(«) + G° s (oj) V sd G d (u) V sd (B3) 



The off-diagonal term G s d describes the interference 
between both channels due to the coupling V s d, and is 
given by: 



G sd {cj) = V sd ■ G° • G d (u>) 



(B4) 



As was shown by Meir and Wingreen in their seminal 
work>Z£, at zero temperature and in linear response the 
exact result for the conductance and current through an 
interacting region reduces to the Landauer result where 
the transmission function can be calculated from the Car- 
oli formula^ 6 . Using the Caroli formula we can calculate 
the coherent transmission via the tip atom as 



T(w) = Tr[T L G\T R G A ] 



(B5) 



where Tl and Tr are the so-called coupling matrices 
describing the coupling of the tip atom A to the left and 
right electrodes: 



with a € {L, R} (B6) 



Hence we find for the transmission T(uj) through the tip 
atom: 




T(w) = T„(u) + Td{co) + T sd {co) 



(B7) 



where T s is the direct transmission through the s-channel 
and Td the corresponding one through the <i-channel 
while T s d describes the transmission involving hopping 
processes between the s- and d-channel: 



T,(u) = 2tt 



T d (w) = 2tt 



(B8) 
(B9) 



T sd (to) - ;l /.I',/,/ + r L>d T R!S ) ■ \G sd \ 2 (BIO) 

In order to calculate the contribution of the s-channel 
to the total transmission, we need to know the spectral 
density of the s-level which is given by the imaginary part 
of the s-level GF given by eq. (|B3[) : 



--Im[G2(w)] • Im[(G°( W )) 2 • G d (u)] 

IT 7T 



(Bll) 



where p° s is the spectral function of the unperturbed s- 
level which in our model is constant, p° s = r s /27r(r^/4 + 
( 2 S ). 5p s (u>) is the change in the spectral density due to 
the coupling to the ti- level with the Kondo peak. For the 
latter we find: 



5p s 



sd 



v. 



{lm[(G^) 2 ] • Re[G d ] + Re[(G°) 2 ] • Jm[G d ]} 



-^■{2- Re[G°] • Im[G°] • Re[G d ] 
(Re[G°] 2 -Im[G°] 2 )-Im[G d ]} 



(B12) 
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The real and imaginary part of Gd are given by: 
z (oj — eji) z x 



ReG d H 
TmG d (w) 



(u-e K ) 2 + T* K T K x 2 + l 

z V k z 1 

{uj - e K ) 2 + T 2 K = IV ' x 2 + 1 



(B13) 
(B14) 



where we have defined x = (ui — ck)/^k- Plugging this 
into eq. (|B12I) we find: 



Sp. 



7T Tk X 2 + I 

■ {2x • ReG° • ImG° - [Re G°] 2 + [ImG°] 2 } (B15) 

We now define the ratio q between real and imaginary 
part of G°: 



ReGg 
ImG" 



2es 



(B16) 



With this we find for the expression in curly brackets in 
cq. (IBTBT) : 

2x • Re G° s • Im G° - [Re G°} 2 + [Im G° s } 2 = (B17) 



[ImG 



0l2 



x t 2x 



ReG° s 
ImG" 



ReG° 



-i 2 



ImG? 



{-2qx-q 2 + l} 



T 2 (1 + q 2 ) 2 

(x + q) 2 - (x 2 + 1) 

r 2 -(i + q 2 ) 2 



Hence we obtain for the change Sp s in the spectral density 
of the s-level due to the coupling to the d-level: 



5p s (to) 



AzV 2 d ( x + q) 2 -(x 2 + l) 



TrT K T2 (x 2 + 1) • (1 + g 2 ) 2 



(B18) 



Summing up, we find for the s-channel contribution to 
the transmission (ll~4"l) : 



T a (u) = T° + 2tt • TL '>_ r *'» . Sps (u) 



1 



2 z V 2 d {x + q) 2 - (x 2 + 1) 
T s -Tk ' (x 2 + l) ■ (l + q 2 ) 



(B19) 



where T" = 2ir p° s ■ (T l.s ■ rfl. s )/r s is the transmission via 
the unperturbed s-channel. 

For the d-channel contribution to the total transmis- 
sion we find according to eqs. (|B7|) and ([3]): 



T d (u) 



• r 



R.d 



1 



r 2 



l 



(B20) 



Finally, for the sd-coupling contribution to the trans- 
mission we need to calculate the absolute square of the 
off-diagonal matrix element G s d of the atomic GF: 



|G sd (u;)| 2 = ^-|G0| 2 -|G d ( W )| 2 
Hence we find: 



v^ 2 



V 2 

1 K 



z 2 V 2 

T sd H=T s °.— ^ 

1 K 



L.d 



R.s 



r 



L.s 



(1+X 2 ) 



1 +X 2 



(B21) 



(B22) 
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